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Abstract 



We study the so-called pir puzzle of the i/j(2S) decay by incorporating 
two inputs; the relative phase between the one-photon and the gluonic decay 
amplitude, and a possible hadronic excess in the inclusive nonelectromagnetic 
decay rate of ip(2S). We look into the possibility that the hadronic excess in 
ip{2S) originates from a decay process of long-distance origin which is absent 
from the J ftp decay. We propose that the amplitude of this additional process 
happens to nearly cancel the short-distance gluonic amplitude in the exclusive 
decay ip{2S) — > 1~0~ and turn the sum dominantly real in contrast to the J ftp 
decay. We present general consequences of this mechanism and survey two 
models which might possibly explain the source of this additional amplitude. 

PACS numbers: 13.25.Gv, 14.40.Gx, 11.30.Hv, 12.39.Mk, 13.40.Hq 
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I. INTRODUCTION 



Absence of the /wr decay mode of ip(2S) has defied a theoretical explanation for more than 
a decade |IJ. The recent measurement by BES Collaboration || has confirmed the absence 
of with even a higher precision, setting its upper bound at a level of factor more than 
60 below what one naively expects from the decay J /if) — » p^ ■ The measurement of other 
decay modes by BES @,[| seems to rule out all possible resolutions for the pn puzzle that 
have so far been proposed by theorists For instance, the large uon branching contradicts 
the helicity suppression |J with or without large intrinsic charm ||. A vector glueball near 
the J /if) mass, if it should exist, can enhance the pir branching for J /if) relative to ip{2S) 
0. However, the magnitude of B(J/if> —>■ pir) is in line with expectation when we compare 
the B{pti)/B{ujti) with the inclusive ratio B(J/if) -> ggg -> X)/B(J/if) -> 7* -> A). What 
happens is not enhancement of pir in J/if> but suppression of pn in ip(2S). [] 

Meanwhile the amplitude analysis of the J/if> decay revealed that the relative phase 
of the gluonic and the one-photon decay amplitude is close to 90° for all two-body decay 
channels so far studied; 1~(T [§, (TO - @, 0, and AA 0. We show in this paper 

that the recent BES measurement in J /if) — > l + _ is also compatible with a large phase. 

In contrast, the pattern of a large relative phase does not emerge for if>(2S). Within 
experimental uncertainties, the relative phase is consistent with zero in the 1~CT decay 
||i2| , [l3H and the l + 0~ decay. This marked difference between J/if) and ip(2S) is another 
puzzle if the three-gluon decay is equally responsible for the strong decay of J/ip and ip(2S) 

There is one more experimental information relevant to the issue. That is the hadronic 
decay rate of if)(2S) which is normally attributed to ip(2S) — > ggg. When we compute with 
the current data the inclusive gluonic decay rate of if){2S) by subtracting the cascade and 
the electromagnetic decay rate from the total rate, it is 60-70% larger, within experimental 
uncertainties, than what we expect from the short-distance gluonic decay alone. This excess 
hadronic branching in if>(2S) may suggest that something more occurs in the gluonic decay 
of ip(2S) than in the J/ip decay. 

In this paper we combine these informations together and search the origin of the marked 
difference between J /if) and if>(2S). While we should be apprehensive about experimental 
errors at present, they might give us a clue to a solution of the pre puzzle. In Section II, 
prompted by the experimental observation in the _ _ , 1 _ 0~, and AA channels, 

we postulate universality of the large relative phase between the gluon and the photon 
decay amplitudes. Specifically, the gluonic decay amplitude acquires a large phase while the 
photon amplitude is real. We point out that a large phase is consistent with new BES data 
in J /if) — > l + 0~. Further progress in the BES analysis in this channel will shed more light. 
We turn to ip{2S) in Section III. The decay branching fractions of ip{2S) — > 1~0~ clearly 
show suppression of the gluon amplitude and favor a small relative phase between the gluon 
and the photon amplitude. We point out that a small phase is more likely in if>(2S) — > l + 0~ 
too. Taking the possible excess in the inclusive hadronic decay rate of ip(2S) seriously, we 
propose that this excess is related to both the suppression and the small relative phase of 
the 1~0~ amplitude. Our proposition is that an additional decay process generating the 



More comparison of experiment with models is found in Ref. 0. 
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excess should largely cancel the short- distance gluon amplitude in the exclusive decay into 
1~CT and that the resulting small residual amplitude is not only real but also destructively 
interferes with the photon amplitude. In Section IV, we first present general consequences of 
the destructive interference. We then examine two scenarios which may possibly generate the 
excess inclusive gluonic decay. One is the contribution of the virtual DD intermediate state. 
The other is a resonance, a glueball or a four-quark resonance, near the ip(2S) mass. Though 
neither idea is novel nor highly appealing, they seem to be among a very few possibilities 
that have not yet been ruled out by experiment. 



II. PHASES OF J/ip DECAY AMPLITUDES 

The relative phase between the gluon and the photon amplitude in the decay J/ijj — > 1~CT 
has been analyzed with broken flavor SU(3) symmetry |§ including the p-u mixing. All 
analyses clearly show that the relative phase should be very large and not far from 90° 
with experimetnal uncertainties. The SU(3) analysis was made also for the 0~0~ modes 
H and the 1~1~ modes JOl for which leading gluon amplitude is SU(3) violating. The 



relative phase were found to be equally large for these modes. Furthermore comparison 
of the electromagnetic form factors in the timelike region with the J/ip decay branching 
fractions revealed that the relative phase is very close to 90° in the NN decay channels too 
|PT| . A question arises as to whether this large relative phase is universal to all decays of 
J/ip or not. There is no persuasive theoretical answer to it at present .0 

In addition to those two-body channels already analyzed, the recent BES measurement 
@ on J/ip -> 1+0", 

B(J/ip -> Kf{im)K T ) = (3.8 ± 0.8 ± 1.2) x 1(T 3 , (1) 
B(J/ip -> Kf{127Q)K T ) < 3.0 x 10~ 3 , 90% C.L. 

is relevant to this issue. We examine here these branching fractions together with the 6 ± 7r T 
branching fraction, B(J/ip -> 6 ± tt t ) = (3.0 ± 0.5) x 10" 3 

Since l£i(1270) and A"i(1400) are superpositions of Ka and Kb of the and the 1#~ 
octet, respectively, 

K+(\m) = K\cos6 + K^smQ, (2) 
^+(1270) = -K\sm6 + K+cos6, 



with 6 m 45° [fL7|, we can parametrize the three decay amplitudes in terms of the gluon 



amplitude a\ of the octet and the photon amplitudes a 7 A/s of an( f 



2 Some attempt was recently made to argue in favor of universal large phases [14|. 

3 A previous analysis [16] of these J/ip decay modes assumed a zero relative phase and used on a 
preliminary value of the upper bound on B(J/tfj — > K^(1270)K T ). Therefore the J/ij) analysis of 
Ref. [16] should be disregarded. However, the analysis of the i^{2S) — * l + _ of Ref. 119] remains 
valid. 
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A(b+ir ) = oi + \/l/5a 7 B 



A(Kf(1270)K~) = (ai + ^l/la lB ) cosO - a 7 ^sin#, 

A(A+(1400)A^) = (ai + ^l/la lB ) sinfl + a^cosfl. (3) 

Since there are two independent helicity amplitudes (or s- and d- waves) for 1 + CT, we should 
use this parametrization separately for the s-wave and the <i-wave amplitudes. The three 
branching fractions can be fitted with 

|ai| > \cl 7 a\ ~ |%b|, 
arg(a^a 7 A) ~ arg(a^a 7B ) « 90°. (4) 

If the ratio T(J/ip — > 7 — ► l + 0~)/T (<//■?/.> — > (7(7(7 — > l + 0~) is comparable in magnitude to 



T(J/tp — > 7 — > X)/F(J/ip — > (7(7(7 — > X) ~ 1/5, we should expect that |a 7j 4/#| ~ 0.7|ai| fTE 
For a 7 ^ = a 7 # = ±0.7mi and ~ 45°, the ratios of the branching fractions prior to the 
phase space corrections, denoted by B , take values as 

fi (& ± 7r T ) : £0(^(1270)1^) = B (K±(U00)K T ) ~ 1 : 0.5 : 0.9. (5) 

While the inequality B{Kf (1270) A^) < B{Kf (1400) A T ) can easily be realized by a wide 
range of parameter values, the other inequality B(b^ 7T^) < 5(^(1400)1^) 

is a little tight. 

If we allow |a 7j 4/B| larger than 0.7|ai| and/or increase the value of 9, however, the current 
central values of the branching fractions are consistent with the large phase hypothesis. We 
should also point out that if the SU(3) breaking correction is made by the meson wavefunc- 
tions C/WU) 2 , it is likely to enhance 5(1^(1400)1^) over S(6 ± 7r =F ) to the direction in favor 
of the large phase fit. 

If we leave a 1 A/B unrestricted in magnitude and phase, a triangular relation holds for 
the amplitudes as 

A(Kf{127G)K*) cos6 + A{Kf (1400) A T ) sinfl = A(bV). (6) 

Determination of the s-to-o? wave ratio of the amplitudes and further study of B(J/ip — > 
K^(14:00)K T ) will eventually resolve the composition of amplitudes and test the large phase 
hypothesis in l + 0~. As it was pointed out previously fll6| , it is also important to resolve the 



discrepancy between B^ir^) and 2B(b 7t ) ||15|| , which theory predicts to be equal. 



To summarize the experimental situation of the two-body J/i/j decay amplitudes, the 
existing data strongly favor large relative phases close to 90° between the gluon and the 
photon decay amplitudes for l - 0~, - - , 1~1~, and NN, and are consistent with a large 
phase for l + _ . 

What does theory say about these relative phases ? In the perturbative picture, the glu- 
onic decay of J/ip proceeds as depicted in Fig. la. The inclusive decay rate is computed with 
the gluons placed on mass shell. In contrast, the photon being far off shell, no corresponding 
on-shell intermediate state appears in the perturbative diagrams of the photon amplitude 
(Fig. lb). Although perturbative QCD (pQCD) is a good description of inclusive charmo- 
nium decays, it is questionable whether it works for two-body decay channels of charmonia. 



To be specific, the pQCD prediction of the asymptotic pion form factor flTJ 
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F^q 2 ) ~ 16na s (q 2 )f 2 /q 2 , 



(7) 



has not been reached at the J/0 mass [19 . Furthermore, the helicity suppression argument 
of pQCD fails for the utt decay channel. Though it is tempting, therefore, we cannot argue 
for the large relative phases on the basis of perturbative diagrams. Whether or not these 
large relative phases are universal to all two-body decay modes of J/ip must be determined by 
experiment. Despite lack of a good theoretical argument at present, we suspect nonetheless 
that the universal large phases so far found are not an accident. 

If the relative phases are close to 90°, it is more likely that the photon decay amplitudes 
are real and consequently the gluon decay amplitudes are imaginary. The reason is as follows: 
In order for the photon decay amplitude to have a substantial phase, the final qq created by 
the virtual 7 should have a large absorptive part of a long-distance origin. This can happen 
if there should be a relatively sharp qq resonance just around J/ip. More likely is that many 
resonances exist below J / ip as in the vector-meson-dominance scenario, or that as in the dual 
resonance model, many increasingly broader resonances appear all the way to high energies 
pO| . In the former case, only the tails of low-lying resonances contribute to the real part. 
In the latter case, a few nearby broad resonances can contribute to the imaginary part, but 
they are outnumbered by many more resonances below and above Jip that contribute to the 
real part. In comparison, we have less insight in hadronization dynamics of the gluon decay. 

Motivated by the results in the amplitude analyses of the two-body of J/tp decays, we 
make two postulates: 

(1) The relative phases between the gluon and the photon decay amplitudes are univer- 
sally large for all two-body decays of J/0. The photon decay amplitudes are predominantly 
real and consequently the gluon decay amplitudes are imaginary. 

(2) The same pattern holds for -0(25) decay as well. 
These are the starting assumptions of our analysis that follows. 



III. rp(2S) DECAY 

A. Relative phase from experiment 

The only large energy scale involved in the three-gluon decay of charmonia is the charm 
quark mass m c . Whether one accepts the argument of the universal large phase in exclusive 
channels or not, therefore, one would naively expect that the corresponding phases should 
not be much different between the J/ip decay and the ip(2S) decay. However, experimental 
data so far available show that the phases are small at least in some two-body decay modes of 
'0(25'). The strongest evidence is in the decay '0(25') — > 1~0~, which includes the puzzling 
0(25) — > p7r. In the case that the final 0~ meson is an octet, we can parametrize the 
0(25) — > 1~0~ decay amplitudes with the SU(3) singlet amplitude ai, the SU(3) breaking 
correction e due to m s — m ud , and the photon amplitude a 1 . The corresponding amplitudes, 
b± and 6 7 , are introduced for the 0~ singlet. For the 0-u; mixing, we assume the nonet scheme. 
The parametrization of the amplitudes [^TJ is listed in Table I for the decay modes so far 



studied in experiment 0. Since the analyzed channels are limited and the uncertainties 
in their branching fractions are still large, we are unable to perform a meaningful x 2 fit at 
present. Therefore, we present only fits to the central values by referring to Table I. 
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First of all, if we ignored the photon amplitude a 7 , we would obtain B(K*°K + c.c.) = 
BiK^K*), which contradicts with experiment, (0.81 ±0.24±0.16) x KT 4 vs < 0.30 x 10" 4 . 
The large splitting between these branching fractions requires that a 7 be comparable to a\. 
If we set a\ and e to zero, we would obtain up to phase space corrections 

B (uon)/B (K* K° + c.c.) = 9/8 (8) 

in contradiction with the measurement, (0.38 ± 0.17 ± 0.11)/(0. 81 ± 0.24 ± 0.16)|] In order 
to come closer to this ratio of the measured values, a large constructive interference should 
occur in K*°K , that is, the relative phase must be small between ai + e and — 2a 7 . Then, 
assuming that ai and e have a common phase, we have a large destructive interference 
between a± and a 7 for both pir and K*^K^ in agreement with experiment.^ This solves the 
pir puzzle and explains also the missing of the K*^K^ mode in experiment. 

With these qualitative observations in mind, we have fitted to the central values of 
the observed branching fractions and then have computed with those parameter values the 
branching fractions of the modes for which only the upper bounds have been determined. 
In Table I we have listed the fit with 5 = — arg(a^a 7 ) = 0° and the large phase fit with 
5 = ±90° for comparison. When the central values of B{wk ) and B(K*°7C + c.c.) are fitted 
with 5 = 0°, the ratio of the photon and the gluon amplitude turns out to be 

a 7 /(ai + e) ~ -0.76. (9) 

For comparison, a 7 /(di + e) ~ 0.14 in the case of J/ip. We would expect |a 7 /ai| ~ 0.22 if the 
ratio of T(ip(2S) -> 7* -> 1-(T) to T(ip(2S) -> ggg -> l~Q-) is roughly equal to T(ip(2S) -> 
7* — > X)/T(ip(2S) — > ggg — > X). Since experiment shows that r(ci;7r)/r(/ + / _ )(oc |a 7 | 2 ) is 
about the same for J ft/) and ■?/>(25'), the large number in Eq.fl^) results from strong sup- 
pression of the total gluonic amplitude a\ + e in ip(2S). As the value of e is varied in the 
range of |e/ai| < 1/3, the value of B(p ± it T ) varies between and 0.04xl0~ 4 . The values for 
_B(p ± 7r =F ) and 5(i^* ± K =F ) can be increased if we stretch within the experimental uncertain- 
ties of B(ujtt) and B(K*°K +c.c). In contrast, the fit with 5 = ±90° overshoots the upper 
bound on B^K^K^) and, if |e| < ||ai|, the upper bound on 5(p ± 7r =F ) very badly. A fit with 
5 = ±90° is virtually impossible even with experimental uncertainties unless |e| ^> |ai|. We 
thus conclude that the relative phase between a\ and — a 7 should be small in ip(2S) — > 1 _ _ 
contrary to the J /if) decay. 

Though it is less conclusive, a small phase seems to be favored in the l + 0~ decay of ip(2S) 
too. It is conspicuous in experiment || that the Kf(lAQti)K T mode is strongly suppressed 
relative to the Kf{l27ti)K^: 

B(iP(2S) -> K±(1270)K T ) = (10.0 ± 1.8 ± 2.1) x 10~ 4 , (10) 
B(ip(2S) -> Kf(U00)K*) < 3.1 x KT 4 . 



It is possible that the SU(3) breaking in the strange and nonstrange meson wavefunctions 
(fnfp/fKfK*) 2 may be responsible for part of the discrepancy. 

5 This main feature of the fit to the i^>(2S) — > 1~0~ amplitudes is found in the earlier paper by 



Chen and Braaten |12|] and, in particular, in the paper by Tuan [13|. 



6 



B(if>(2S) — > ^ti^) is half way between them ||: 

B(i/j(2S) -> fcV) = (5.2 ± 0.8 ± 1.0) x 10~ 4 . (11) 

We can use Eq. (4) as the parametrization of if)(2S) — > l + 0~. First of all, if a\ dominated 
over CL^A/Bi we would have B{bV) ~ 25(7^(1270)/^) ~ 25(^(1400)^) for ~ 45° 
in disagreement with experiment. Just as in ip{2S) — > 1~0~, |ai| is comparable to |o 7j 4/b|- 
Next, the strong suppression of K^(14:00)K T relative to K^(1270)K T can be realized only 



when ai + yl/5a 7 s interferes destructively with a^A- Therefore, the relative phase between 
di and a^A/B must be small modulo it. To obtain B[Kf (1270) K T ) « 2i?(6 ± 7r T ), we need 



a\ + y l/5a 7 B ~ — a 7 A- The allowed range of the amplitude ratios was plotted in Ref. [TJ 
by choosing all amplitudes as relatively real and assuming tentatively the s-wave decay for 
phase-space corrections. Though it is not impossible to fit the three branching fractions 



with 9 ~ 90°, we must have a\ ~ and y l/5a 7B ~ — a 7 A in that case. 

To summarize for ip(2S), the data on 1~0~ virtually excludes the possibility of a small 
phase between a% and a 7 . A fit to l + _ has more room when the relative phase is small. 
There is no evidence for that the relative phase must be large in ip(2S). For some reason, 
a large relative phase does not seem to occur in the ip(2S) decay. We ask what causes this 
marked difference between J ftp and ip{2S) when we postulate the universal large phase for 
ip{2S) as well as for J /if). 

B. Excess hadronic rate in inclusive hadronic decay 

It has been noticed that when one computes the inclusive hadronic decay rate of if)(2S) 
through ggg by subtracting the rates of the cascade and electromagnetic decays from the 
total decay rate, it is substantially larger than what we expect from an extrapolation of J /if). 
The number with a conservative error estimate based on the listings of Reviews of Particle 
Physics |T5 



is 



B(ib(2S) -> ggg + 007) , s 

X, TI — yyyjLMiiL = 0.23 ± 0.07, 12 
B(J ip^ggg + ggj) 1 ; 



which should be compared with 

'a s {if){2S))\ 3 B{tP{2S) ^l+i 



a s (J/if)) J B(J/if) 1+1-) 



= 0.134 ±0.034. (13) 



Smaller errors (0.226 ± 0.052 vs 0.141 ± 0.012) have been attached in a recent literature |2B| 
with a different error estimate. We would expect that the two numbers should be equal to 
each other since the wavefunctions at origin appear in common in Eqs. (|T2"D and (13). The 



discrepancy of 60-70% between them alarms us particularly because all numbers involved 
have been repeatedly measured over many years.f] 



6 The author learned that BES collaboration is considering a different determination of the cascade 



decay branchings [24] 



7 



In comparison we find no similar excess in T(2S) though experimetnal undertaintues 
are large. In terms of the ratio of branching ratios, B(ggg + gg^f) = B(T — > ggg + 
gg7)/a s (Y) 3 I?(Y — > three T's are more in line: 

_ f (4.5 ±0.2) x 10 3 , T(IS) 

B(ggg + ggi) = (4.9 ± 0.9) x 10 3 , T(2S) , (14) 
[ (4.0 ±0.4) x 10 3 , T(3S). 

where the total leptonic branching for T(3S)Q has been substituted with three times 
B(fi + fi~), the only quoted leptonic branching. It appears that the excess in B(ggg + ggj) 
is unique to tp(2S). However, this excess in the inclusive rate has not shown up in the rates 
of the exclusive channels so far measured. In fact, the ratio B(ip(2S) — > h)/B(J/tp — > h) 
scatter around the expected value (~ 13-14% of Eq.(|i~3"D), which was often called the 14% 
rule. Some remarks should be in order on it. 

First of all, the 14% rule is largely violated in many of two-body and quasi-two-body 
channels, as we recently learned in the BES data 0. The pir channel is an extreme case. For 
multihadron channels, there are actually not so many modes that are available for testing 
the 14% rule.Q In Table II, we have tabulated the ratios for the modes not listed in 
but available for comparison. We see that the ratios scatter rather widely above and below 
14% with some tendency of being smaller than 14%, but with fairly large experimental 
uncertainties. It is important to notice that the branching fractions of all modes in Table II 
add up to no more than 15% of the total gluonic decay branching of J /if). Indeed, only one 
charge state has been available for comparison from each of 5n, 7n, 2nKK, and NNnn. We 
have not yet seen comparison of the rest. The so-called 14% rule is based on very limited 
number of decay modes. It is premature to preclude the hadronic excess with the data of 
multihadron exclusive channels. 

If future experiment shows that the hadronic excess in ip{2S) is real, it may have some- 



thing to do with the /wr puzzle p3] and with the abrupt change of the relative phase of 
amplitudes from J/ip to ip(2S). A process responsible for the excess inclusive hadron rate 
can interfere with the short-distance gluon decay in exclusive modes. If its amplitude makes 
a large destructive interference with the three-gluon amplitude in ip(2S) — > 1~0~ and if the 
sum is nearly real and comparable to the photon amplitude in magnitude, our puzzle can 
be solved. We shall look into possible sources of this rate excess in the following. 



The branching to T (25)77 quoted in [15] is sum of the cascade decay branchings Xbj{^P)l - > 
T(2S')7 in view of the 77 invariant mass spectrum |25j and also of its magnitude (> S(T(25)7r°7r )). 



It is counted in the radiative decay branchings separately listed in 



Gu and Li lumped multihadron modes together and compared between J/ip and tp(2S). Then 
the number is dominated by three modes, 7r + 7r _ 7r°, 2(7r + 7r _ )-7r and 3(7r + 7r~)7r°, which happen to 
be below 14%. The mode tt + tt^tt° is actually pix and its nonresonant content is consistent with 



zero [15]. More revealing are the ratios of individual modes. 



S 



IV. ADDITIONAL HADRONIC AMPLITUDE IN ip(2S) 



A. General consequences 

If the origin of the problem is in the interference of an unknown additional process with 
the short-distance gluon decay of ip(2S), we expect a general pattern of correlation between 
the decay angular distribution and suppression or enhancement. 

The decay angular distribution for two final hadrons is generally of the form, 

dT/dtt oc 1 + a cos 2 9, (\a\ < 1) (15) 

where 9 is the polar angle measured from the e + e~ beam direction. For 1~0~ and 0~0~ 
decays, the value of a is constrained kinematically to +1 and —1, respectively, while it is 
determined dynamically by the helicity content, ±1 or 0, of the final state in other decays. 
In 1~0~ and 0~0~, therefore, any additional amplitude has the same angular dependence 
as the three-gluon and the photon amplitude irrespective of its origin. Consequently a high 
degree of destructive or constructive interference with an additional amplitude is possible 
in these decays. Observation of the strongest suppression in the 1~0~ mode is consistent 
with this pattern ||. We expect that the decay rates of ip{2S) — > CTCT may also be quite 
different from those of J /if) -> (TO - . In terms of B(0 _ _ ) = B(Q-Q-)/oP 8 B{n + n~), the 



current data jnj give 



■n, + S I 0.15 ±0.02 for J/4 

S(7F 71 H 0.8 ±0.5 for i)(2S) (16) 



B(K + R- 



0.24 ± 0.03 for J/ijj 
0.94 ± 0.66 for ijj(2S) 



Within the large experimental uncertainties we see a hint of large constructive interference in 
ip{2S) — > 0~0~. In contrast, in other processes an additional amplitude and the three-gluon 
amplitude have different angular distributions in general. A large interference can occur only 
when dynamical mechanisms of two processes are similar. Otherwise it should be a result 
of a high degree of accident. When a large disparity is observed between corresponding 
two-meson decay rates of J/ip and ip(2S), therefore, the decay angular distribution of this 
channel will also be very different between J/ip and ip(2S). This will give a good test of the 
idea of interference with an additional amplitude. 

The other consequence is in multibody final states. Since the additional process enhances 
the inclusive rate, a large number of exclusive decay channels should receive enhancement 
rather than suppression. When there are many hadrons in the final state, chance of interfer- 
ence between amplitudes of different decay mechanism is much smaller because of difference 
in subenergy dependence and event topology. Therefore enhancement will not be dramatic. 
While we have not yet seen such enhancemnet in Table II, we expect that the branching 
fraction tend to be enhanced in many nonresonant multibody channels of ip(2S) relative to 
J/v. 

Where does the additional amplitude possibly come from ? There are few options left in 
modifying charmonium physics radically. Since aspects of perturbative QCD have been well 
understood, we are bound to look for the origin of the problem in long-distance physics of 
one kind or another. 
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B. ip(2S) -> DD -► hadrons 

One unique feature of if){2S) is a close proximity of its mass to the -D.D threshold. The 
^(25) mass is only 43 MeV (53 MeV) below D°D° (D+D~), while T(3S) is 200 MeV 
away from the BB threshold. Can the small energy difference^ between if)(2S) and DD 
have anything to do with the excess ? It may happen that if)(2S) picks up a light quark 
pair through soft gluons and dissociates virtually into DD, which in turn annihilate into 
light hadrons. (See Fig.2.) The dominant process of the DD annihilation is through cc 
annihilation through a single hard gluon. The small energy denominator enhances creation 
of virtual DD while p-wave creation compensates the enhancement. It is difficult, actually 
nearly impossible, to give a reliable computation of this sequence. Deferring estimate of 
the rate to future, we shall comment here only on whether the DD contribution can have 
a final-state interaction phase large enough to cancel the perturbative gluon amplitude or 
not. 

Since the DD intermediate state is above the if>(2S) mass, a phase of amplitude must 
come from the subsequent annihilation of cc and thereafter. After an energetic light quark 
pair qq is created from cc, each of qq picks up a soft light quark from the light quark cloud of 
DD to form mesons. (See Fig.3a.) Kinematically, this step of the hadron formation process 
is quite different from that of the timelike electromagnetic form factor of a meson in which 
energetic light quarks pick up collinear quarks created by a hard gluon. (See Fig. 3b.) In our 
case color-dipole moment is large for all pairs of quarks p7 |. Furthermore, the cm. energy 



of a hard quark in one meson and a soft quark in the other meson is in the low energy 
resonance region, 



= 0{pK QCD m c ) < lGeV. (17) 

Therefore, one cannot argue that final-state interactions should be small between final 
mesons. We expect that there is a good chance for the amplitude of ip{2S) — > DD — > mesons 
to acquire a substantial final-state interaction phase. Weakness in this argument is that the 
phase can be large but need not be large. 

The idea of the virtual DD dissociation actually has some common feature with that of 
the "higher Fock component" of charmonia [p8p9f. The DD state can be viewed as part 
of the four-quark Fock space of ip(2S). The higher Fock component was proposed as an 
additional contribution to J/if) — > 1~0~ to solve the pn puzzle j29[| . It was argued that it is 



more significant in J/if> than in ip(2S). As we have emphasized, however, there is nothing 
anomalous about J /if) ~ > 1~0~. 

C. i/)(2S) -^resonance^ hadrons 

The second idea is a twist of an old one: A noncharm resonance may exist near the 
if){2S) mass and give an extra contribution to the hadronic decay rate. A glueball was 
proposed earlier at the J/ip mass to boost the decay rate of J /if) M. However, we now 



9 This was brought to the author's attention by J.L. Rosner [26| 
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want it near ip{2S) not near J/ip- We look into the possibility that some resonance around 
the if}(2S) mass, a glueball or four-quark, destructively interferes with the perturbative 
ip{2S) — > ggg — > 1~0~ decay. Admittedly, the idea is ad /ioc and there is some difficulty 
aside from unnaturalness. 

Light-quark resonances of high mass (~ 3.7 GeV) and low spin are normally too broad 
to be even recognized as resonances. Four-quark resonances may be an alternative if they 
exist at all. The mass of 3.7 GeV is normally considered as too high for the lowest vector 
glueball. An excited glueball state of J P C = 1 serves our purpose. Whatever its origin 
is, let us introduce here such a resonance, call it R, and see its consequences. 

In order for ip(2S) to decay through R as strongly as through three gluons, the coupling 
/ of R to if)(2S) defined by —fiJiRip^R^ must be large enough. The if)(2S)-R mixing at the 
ip{2S) mass is given by 

1 (18) 



Am — iT 



R 



where Am = 771r — m(ifj(2S)) and is the total width of R. It leads to T(tfj(2S) — ► 
R — > hadrons) ~ \f\ 2 /Tn when |Am| < 0(Tr). To obtain T(ip(2S) — > R — > hadrons) ~ 
T(ip(2S) — > ggg), we need therefore 

\f\ 2 ^T R T^(2S)^ggg). (19) 

If R is a light-quark resonance qq, \ f\ would be much too small for the following reason: While 
\f\ 2 is of the order of T(ip(2S) — > ggg)T(R — > ggg) for qq, we expect T(R — > ggg) <C T R 
because of the suppression of qq — > ggg. Therefore there is no chance to satisfy Eq. (|T9|) . 
It is likely that the same argument applies to four-quark resonances. 

For glueballs, we simply do not have enough quantitative understanding to rule out a 
large enough coupling to ip{2S). Hou and Soni [[?] proposed a glueball near the J/ip mass in 
order to enhance J/ip pn (and its symmetry- related modes) but not other decay modes. 
To accomplish it, this glueball must have very special, if not unnatural, properties |30|: It is 



nearly degenerate with J/ip with a quite narrow width for an object of mass ~ 3 GeV and 
decays predominantly into 1~0~. Later Hou [[J!]] relaxed the constraint on the p7r branching 
to argue that such a glueball was not yet ruled out by the search of the BES Collaboration 



In our case, since a glueball is introduced to account for the hadronic excess, it should 
couple not primarily to the 1~0~ channels, but to many other channels. What we need is 
a generic vector glueball with no special or unusual properties. If V R is as narrow as 100 
MeV, for instance, the mixing |e| = O(10 -2 ) would be able to account for the excess in the 
inclusive hadron decay of ip(2S). The width can be wider. In that case the pole transition 
strength / should be stronger according to Eq . fll9D . Since the glueball R couples to a photon 
only indirectly through its mixing to a quark pair, it is hard to detect R in the hadronic 
cross section of e + e~ annihilation near the ip(2S) mass. Searching by hadronic reactions 
such as pp annihilation is a daunting task. From a purely experimental viewpoint, such a 
resonance has not been ruled out |Hfl . 



Though the resonance scenario is admittedly a long shot, it is one of a very few options 
left to us. One reason to pursue this somewhat unnatural scenario is that the amplitude for 
ip{2S) — > R —>■ X has automatically a large phase when Am < 0(T R ) since the coupling 
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/ is likely real, that is, dominated by the dispersive part. If that is the case, the resonant 
amplitude can interfere strongly with the three-gluon amplitude in two-meson decays. 

V. CONCLUSION 

We have searched for a clue to solve the pn puzzle in this paper. Our purpose is to locate 
the source of the problem rather than to offer a final solution of the problem. Two threads 
have been exposed which may eventually lead us to a solution of the pn puzzle. They are 
the phases of the decay amplitudes and a possible excess in the inclusive hadronic decay rate 
of ijj(2S). An experimental confirmation of the excess will be the most useful in directing 
theorists. If it is confirmed, it will be quite an important experimental discovery by itself. 
One crucial experimental information will be the angular distributions of J/ip and ip(2S) 
into the channels other than 1~0~ and 0~0~. Difference in the angular distributions should 
have direct correlation with enhancement and suppression in general. As for the source of an 
additional process, the virtual DD pair and the vector glueball are two options that cannot 
be ruled out. To be frank, however, we admit that both ideas have unnaturalness. More 
an attractive alternative is highly desired. It is possible that the large relative phases so far 
observed in J/ip decay are an accident and that the pn puzzle is a problem of incalculable 
long-distance complications. However, our hope is that there might be something novel, 
simple, or fundamental hidden beneath the issue. 
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TABLES 



TABLE I. Parametrization of the ip(2S) — > 1~0~ amplitudes and values of branching fractions. 
The amplitude e represents the T33 breaking of m s — m U( i instead of As breaking. The 77-7/ mixing 
angle is chosen to be Op = —20°. The fits to the central values are shown for the minimum and 
the maximum relative phase, 5 = — arg(aja 7 )= 0° and ±90°. The ranges of values for pir and lot] 
are given for \e/a\\ < 1/3. 



Modes 


Amplitudes 


Branchings (in 10 4 ) 


Fits 


6 = 


5 = ±90° 


p+Vr-(= p°7T ) 


a\ + a 7 


< 0.09 


0~0.04 




0.21~0.70 


K*+K- 


ai + e + a 7 


< 0.15 


0.00 




0.30 


K*°K° 


ai + e — 2a 7 


0.41 ±0.12 ±0.08 


0.41 




0.41 


LOTT 


3a 7 


0.38 ±0.17 ±0.11 


0.38 




0.38 


UJT] 


V /T7 3 "( a i + a ~f) ( W7 ?8) 


< 0.33 


0.06 ~ 0.22 




0.05 ~ 0.31 


uorj 


V / 273(6i + fe 7 ) {corn) 


0.76 ±0.44 ±0.18 


0.76 




0.76 



TABLE II. Branching fractions of J/ip and ip(2S), and the ratio B(ip(2S) -> h)/B(J/ip -> h). 
The 7T + 7r~7r mode is not included here since J/ip — ► /t is entirely J/ip ^ pit within experimental 
uncertainties. 



Modes 


Branchings of J/ip 


Branchings of ip{2S) 


Ratio 


2(7r + vr-7r°) 


3.37 ±0.26 x 10~ 2 


3.0 ±0.8 x 10" 3 


8.9 ±2.5 % 


3(7r + 7r-)7T° 


2.9 ±0.6 x 10~ 2 


3.5 ± 1.6 x 10 -3 


12.1 ±6.1 % 


K+K- 


2.37 ±0.31 x 10~ 4 


1.0 ±0.7 x 10~ 4 


42 ± 30 % 


■K+iT-K+K- 


7.2 ±2.3 x 10~ 3 


1.6 ±0.4 x 10~ 3 


22.2 ± 8.8 % 


pp 


2.12 ±0.10 x 10 -3 


1.9 ±0.5 x 10" 4 


8.9 ± 2.4 % 


ppir 


1.09 ±0.09 x 10 -3 


1.4 ±0.5 x 10~ 4 


12.8 ± 4.7 % 


PP1T + TT~ 


6.0 ±0.5 x 10~ 3 


8± 2 x 10~ 4 


13.3 ±2.1 % 
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FIGURES 




FIG. 1. Decays of charmonium into two mesons (a) through ggg and (b) through one pho- 
ton. The vertical broken line in diagram (a) indicates that the gluons are placed on the mass 
shell when the inclusive decay rate is computed with ggg.li perturbative QCD dominated, the 
one-gluon-exchange diagram depicted in (b) would dominate in the final state of the one-photon 
process. 




FIG. 2. The decay ip(2S) — > DD (off shell) — ► meson + meson. 
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FIG. 3. (a) Formation of a light meson pair by energetic qq from cc and wee qq in DD annihi- 
lation. The invariant mass is small for qq and for qq. The arrows denote directions and magnitudes 
of momenta, (b) Light meson pair formation in one-photon annihilation where all quarks are hard. 
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